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The de Gennes’ blob model is extensively used in different problems of polymer physics. This 
model is theoretically applicable when the number of monomers inside each blob is large enough. 

For confined flexible polymers, this requires the confining geometry to be much larger than the 
monomer size. In this manuscript, the opposite limit of polymer in nano-slits with one to several 
monomers width is studied, using molecular dynamics simulations. Extension of the polymer inside 
nano-slits, confinement force on the plates, and the effective spring constant of the confined polymer 
are investigated. Despite of the theoretical limitations of the blob model, the simulation results are 
explained with the blob model very well. The agreement is observed for the static properties and 
the dynamic spring constant of the polymer. A theoretical description of the conditions under which 
the dynamic spring constant of the polymer is independent of the small number of monomers inside 
blobs is given. Our results on the limit of applicability of the blob model can be useful in the design 
of nano-technology devices. 


I. INTRODUCTION 


Advances in nanotechnology have enabled investiga¬ 
tions of single polymers confined in the nano-scale [l|,I2- 
This problem has attracted interests due to its applica- 
tions in the design of nanotechnology devices [^-Q- It 
also helps in obtaining more knowledge of polymers con¬ 
fined in biological environments 0- Besides, it has 
become possible to check accuracy of the old theories 
and improve them Ir |9||. Single polymers confined in- 
side nano-spheres [ifl. 1llj|. nano-channels 0, f 
and nano-slits (between two parallel plates) [1 
are studied. 

The most known theory for confined polymers inside 
nano-channels or nano-slits is the de Gennes’ blob model. 
For semi-flexible polymers, Odijk theory is relevant when 
the size of the confining geometry, D, is smaller than the 
persistence length of the polymer, Lp i[ll- Recently, it 
is shown that extended de Gennes’ theory should be used 

in the case of semi-flexible polymers when Lp < D < — 
0. Here, w is the monomer width. 

Most of the studies on confined polymers focus on 
double-stranded DNA, which is a semi-flexible polymer. 
Recently, single stranded DNA is introduced as a model 
system for study on flexible polymers [^ . Also, many 
proteins and synthetic polymers are flexible. For flexible 
polymers, the only relevant theory is the de Gennes’ blob 
model and it would be useful to better understand its 
limitations and regime of applicability. The blob model 
has broad applications in different problems of polymer 
physics, such as polyelectrolytes [2l|, polymers under 
tension and polymers adsorbed to surfaces or 
confined near surfaces unfer external fields [2^. 

In the blob model, the confined polymer is divided into 
smaller sections, called blobs, that have the size of the 
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confining geometry, D. Inside blobs, polymer does not 
feel that it is confined. The number of monomers inside 

each blob is found from the relations of the self-avoiding 

2 _ 

walk (SAW) in free space, g ^ (t)"- H^re, v is the 
Flory exponent in 3D and b is the monomer size. In 
length scales larger than the blob size, the polymer (with 
N real monomers) behaves like an effective polymer with 
y monomers of size D. This effective polymer has a 
SAW in one or two dimensions for a nano-channel or a 
nano-slit, respectively. The number of monomers inside 
each blob, g, and total number of the blobs, should 
be large enough, to be able to use the relations of sta¬ 
tistical mechanics for SAW inside or outside the blobs, 
respectively (^ . 

A number of studies on confined polymers are per¬ 
formed in conditions that the above two requirements are 
not satisfied. It is shown that the blob model is applica¬ 
ble in micro-devices, where the number of blobs is very 
small [1^. Also in some practical situations, the poly¬ 
mer is confined in severe conditions that the number of 
monomers inside each blob is very small . Specially in 
nano-devices for DNA barcoding, full stretching of DNA 
is required 0 . Computer simulations show that the blob 
model can explain the statics of a flexible polymer inside 
narrow channels [HQ. However, the polymer dynam¬ 
ics in nano-channels is only described with the blob model 
when the number of monomers inside each blob is large 
enough MM- 

In this manuscript, a flexible polymer confined inside 
a nano-slit is studied using molecular dynamics (MD) 
simulations. The extreme condition that the number of 
monomers inside each blob is small is considered. It is 
shown that our results for the radius of gyration and con¬ 
finement free energy of the polymer can be explained very 
well, using the blob model. Despite previous studies (Ref. 
[l^ and references therein), no correction is made to the 
slit width (distance between the two plates). Agreement 
between simulation results and the blob model is ob- 






2 


served even for very narrow slits with 1-2 nionomer(s) 
width. These widths are completely below the regime of 
applicability of the blob model. 

Extension of the polymer inside the nano-slit has fluc¬ 
tuations around its mean value. From these fluctuations, 
an effective spring constant can be obtained. It is a dy¬ 
namic variable for the polymer and is a measure of elas¬ 
ticity of the polymer. Our results for the effective spring 
constant are in good agreement with the blob model, un¬ 
like the case of the polymer inside the nano-channel [H- 

M- 

This discrepancy is due to the difference between the 
origins of the effective spring constant of a polymer con¬ 
fined in a nano-channel and a nano-slit. Two factors 
govern the fluctuations of a polymer confined in a nano¬ 
slit: fluctuations in the size of blobs and fluctuations in 
the self-avoiding walk of the blobs. It is obvious that the 
latter factor is dominant. However, for a polymer con¬ 
fined in a nano-channel, the blobs are arranged one after 
the other along the channel. So, fluctuations in the blob 
size are the only source of fluctuations in the extension 
of the polymer (see Fig. [I}. This makes the internal 
dynamics of the blobs important, in the dynamics of the 
whole polymer. As a result, it is determining to have 
the sufficient number of monomers inside each blob, for 
a polymer confined inside a nano-channel. 

In the next section, the statics and dynamics of a flex¬ 
ible polymer confined in a nano-slit is reviewed briefly. 
In Sec uni our simulation method is described. Results 
of the simulation and their compatibility with the the¬ 
ory are discussed in Sec. IIVI Finally, a summary of the 
results is presented in the last section. 


II. THEORY 


Consider a flexible polymer confined between two par¬ 
allel plates (Fig. [IJ. Confinement extends the polymer 
on the walls, to a radius of gyration larger than its free 
radius of gyration Rg ^ bN‘^. Here, N is the number 
of monomers of the polymer, b is the monomer size and 
iz = I is the Flory exponent in 3D. Sections of the chain 
smaller than the distance between the plates, D, do not 
feel the constraints. These sections, called blobs, have a 
statistics similar to a chain in free space D ^ bg's, where 
g is the number of monomers inside each blob. 

The polymer is composed of ^ blobs that cannot pen¬ 
etrate into each other. So, it can be regarded as an effec¬ 
tive two-dimensional polymer with y monomers of size 
D. We use the exponent of a self-avoiding walk in two 
dimensions, v = j, to find the extension of the polymer 
parallel to the plates i?|| D{^)i. This gives the rela¬ 

tion between the average extension of the polymer with 
the total number of monomers and the distance between 
the plates, 




(a) 



FIG. 1. (Color online) (a) Blob model for a flexible polymer 
confined inside a nano-channel. The blobs of the polymer in 
the nano-channel are arranged one after the other along the 
nano-channel, (b) Top view of a polymer confined inside a 
nano-slit. The blobs of the polymer in the nano-slit span the 
2D space between the plates. The number of monomers inside 
one blob and total number of the blobs are clear in the two 
cases (a) and (b). 


i?|j ^ bN^ 



( 1 ) 


Free energy of confining the polymer is ksT per blob, 
Fconf ~ ksT^. Substituting g in this relation gives 


Fconf 


ksTN 



( 2 ) 


The force exerted on the plates is found by differenti¬ 
ating the confinement free energy with respect to D, 


/ 



( 3 ) 


The polymer has fluctuations around its average ex¬ 
tension. The mean square deviation from the average 
extension is described by an effective spring constant 
(Ai?^|) = {{R\\{t) - R\\avY) ~ ksT/keff. The simplest 
way to calculate the spring constant, kefj, is to elimi¬ 
nate D between the relations [T] and [51 This gives the free 
energy as a function of extension. Differentiating the 
polymer free energy twice with respect to its extension 
and substituting the equilibrium extension gives a spring 
constant Q. This calculation does not give a correct re¬ 
sult, because equation is the free energy of the polymer 
at its equilibrium extension. Indeed, the free energy of 
the polymer for any desired extension should be used in 
the calculation. 

Suppose the two ends of the confined polymer are 
under tension and it is elongated to the extension Rf. 
Free energy of a two- or three-dimensional polymer with 
N monomers of size b which is extended to size Rf is 
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Ftens ~ ksT [l^- V is equal to | for 3D and 

equal to I for 2D. By differentiating this free energy twice 
with respect to Rf and substituting the equilibrium ex¬ 
tension Req ^ bN'', the spring constant is obtained as 
keff = ~ ~ This is in agree¬ 

ment with the relation used in Ref. 1^. 


Here, the confined polymer is described as an effec¬ 
tive polymer with monomers of size D in two di¬ 
mensions {v = |). So, the spring constant becomes 
keff ~ 7 — ^^*^ 3 ' 2 ■ Substituting for g, we get the spring 

constant for the confined polymer 


62 ^ U; ^11 ■ 


(4) 


In the final relation, i?|| is substituted from eq. [T] To 
derive this relation, the confined polymer is simply mod¬ 
eled as an effective two-dimensional polymer composed of 
blobs and the effect of the external tension is to move the 
confinement blobs toward a linear alignment. This simple 
description is valid for small tensions with Rf < D(y). 

The blob model is accurate when the number of 
monomers inside each blob, g, and total number of the 
blobs, y, are considerable. This is necessary to be able 
to use statistical mechanics relations for a self-avoiding 
walk inside and outside the blobs, respectively. The size 
spanned by a self-avoiding walker in space is proportional 
to the number of steps to the power of the Flory exponent 
in two or three dimensions, when the number of steps is 
large enough. For g to be large enough, the confining ge¬ 
ometry (nano-slit or nano-channel) should be much larger 
than one monomer size. 

It is instructive to compare the effect of an external 
tension on a polymer confined inside a nano-channel or a 
nano-slit. Prior to tension, the polymer inside the nano¬ 
slit arranges its blobs in two dimensions. The effect of the 
external tension is to move the blobs toward a one dimen¬ 
sional arrangement. However, the blobs of the polymer 
inside a nano-channel have already lost two degrees of 
freedom and lie on a one dimensional line. So, the poly¬ 
mer has to stretch its blobs to respond to the external 
tension. One can conclude that the elasticity of a poly¬ 
mer inside a nano-slit has its origin in the arrangement 
of the blobs. This is in contrast to the elasticity of a 
polymer inside a nano-channel which originates from the 
elasticity of each blob. As a result, to use the blob model 
for the former case, the number of blobs should be suf¬ 
ficient; but for the latter case, the number of monomers 
inside each blob is also important. 


III. SIMULATION METHOD 

We use MD simulations to check the limits of accu¬ 
racy of the blob model for a flexible polymer confined 



FIG. 2. (Color online) Side view of a polymer confined in a 
nano-slit. The surface of each plate is defined such that the 
closest distance of the center of a monomer to the plate is 
equal to 6/2. Distance between the plates in our paper, D, 
is shown in the hgure. An effective width for the nano-slit is 
also shown, to examine the effect of adding correction to the 
slit width. 


in a nano-slit. A schematic of the system is shown in 
Fig. [21 The polymer is modeled as a bead-spring chain. 
Monomers of the polymer are connected by the FENE 
potential. 


UpENEir) 


-^Krl In 
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(5) 


Shifted-truncated Lennard-Jones is used for the excluded 
volume interactions between the monomers and between 
the monomers and the plates. 


ULj{r) = 4e 




T <2^ b. (6) 


e and b are the energy and length scales of the simulation. 
K = lOOe and tq = 1.56 are used for the FENE potential. 

The plates are infinite and the closest distance of the 
center of a monomer to the plate surface is 6/2. This 
means that the spherical surface of the monomer can 
touch the plate, in agreement with our intuition from 
a wall (Fig. |2]). In previous studies, the center of a 
monomers was allowed to touch the impenetrable plate 
surface [l^. So, it was needed to add the monomer size 6 
to the distance between the walls to compare the results 
with the blob model. It is important to note that this 
is not attributed to the finite size of the chains in these 
simulations. Here, it will be shown that the simulation 
results for similar chains are in agreement with the blob 
model, without any corrections in the wall distance. 

The equations of motion are integrated using the Ve¬ 
locity Verlet algorithm, with the step size equal to O.OItq. 

To = is the MD time scale and m is the monomer 

mass. The system is kept at the constant temperature 
T = using the Langevin thermostat with the fric¬ 

tion coefficient I. OtT ^. The simulations are performed 
using ESPResSo [^ . 
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FIG. 3. (Color online) Simulation results for i?|| versus the 
theoretical relation J\f0.75j-^-o.25^ _l_^ ^ 

and l> show the simulation results for D/b=l-9, in ascending 
order. Inset: Distribution of the monomers in the distance 
between the plates. The solid, dashed and dash-dotted lines 
correspond to the three narrowest slits D/b=l-3, in ascending 
order. 


In our simulations, the nano-slit width is changed from 
b to 96. Polymers with 100, 200 and 300 monomers are 
investigated. For slits with width 3-96, the monomers are 
arranged in an initial structure close to a self-avoiding 
walk. Then, the system is warmed up to remove any 
overlap between the monomers. For slits with width 6 
and 26, the monomers are initially arranged on one line. 
This initial configuration reaches equilibrium in a longer 
time interval, but it is inevitable because of the very small 
width of the slits. 

All simulations are performed for a total time of 
Data points to find time averages of the static and dy¬ 
namic properties of the polymer are collected after the 
time in time intervals equal to tq. In all plots, the 

size of errorbar at each point is smaller than the symbol 
size. 

The radius of gyration of the polymer parallel to 
the plates is calculated using the relation i?j| = 

XcmY + (Yi - Ycm)^- Xi and Xcm are the 
X components of the positions of the ith monomer and 
the center of mass, respectively. Yi and Ycm are the 
corresponding y components. The force on each plate in 
all steps of the simulation is sum of the Lennard-Jones 
forces from all monomers interacting with the plate. 



log (N) 

FIG. 4. (Golor online) logio-logio plot of radius of gyration 
of the polymer parallel to the plates i?|| versus the polymer 
length, N, for different nano-slit widths. O, □, x, •, -I-, 
V, A, <1 and [> show the simulation results for D/b=l-9, in 
ascending order. The corresponding a values are 0.74, 0.78, 
0.79, 0.76, 0.72, 0.74, 0.74, 0.78 and 0.79. The average of a is 
0.76. 


IV. SIMULATION RESULTS 
A. Radius of gyration 

Radius of gyration of the polymer parallel to the plates 
versus the theoretical relation in equation [T] is shown in 
Fig. [3] All data lie on a single line, in the whole range 
of parameters studied. This shows agreement between 
the simulation results and the blob model. Distribution 
of the monomers in the distance between the plates is 
shown in the inset of Fig. [31 for the three narrowest slits. 
Because of the soft nature of the potential, monomers 
can penetrate the walls rarely. 

To further investigate the correspondence of the sim¬ 
ulation results with the theory, exponents a and j3 in 
the relation i?|| oc N°‘D^ are obtained from simulation. 
These exponents are found from the slopes of the logio- 
logio plots of the radius of gyration versus the polymer 
length and the nano-slit width, respectively. As can be 
seen in Fig. |4l a is close to the theoretical value even 
for very narrow slits (1-26). a is not a decreasing or in¬ 
creasing function of the nano-slit width. So, the small 
number of monomers inside each blob is not determining 
in this exponent. The average of a is 0.76, which has %1 
deviation from the theoretical value. Fig. EKa) shows the 
exponent /3 for different lengths of the polymer. It is seen 
that finite size of the chains has no effect on the simula¬ 
tion results. The results are summarized in the relation 
i?ll oc iV0'76i:)-0-25. 

For comparison with Ref. [T^ . an effective width is 
defined for the nano-slit by adding a monomer size to the 
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log [(D+5)/b] 

FIG. 5. (Color online) logio-logio plot of i?|| versus the nano¬ 
slit width, (a) and (b) show two different definitions for the 
nano-slit width (see Fig. n. As can be seen, omitting 5 results 
in an excellent improvement in the scaling exponent. 


nano-slit width (see Fig. [2|) in Fig. [SKb). The exponent 
P deviates largely from the theoretical value, with this 
change in the dehnition of the wall surface. 



FIG. 6 . (Color online) Simulation results for / versus the 
theoretical relation ND~^'^. O, □, x, •, -f, V, A, <1 and > 
show the simulation results for D/b=l-9, in ascending order. 



log ( N) 


FIG. 7. (Color online) logio-logio plot of / versus the polymer 
length, N, for different nano-slit widths. O, □, x, •, -I-, 
V, A, < and > show the simulation results for D/b=l-9, in 
ascending order. The corresponding 7 values are 1, 0.99, 0.99, 
0.98, 0.99, 0.96, 0.98, 0.95 and 0.99. The average of 7 is 0.99. 


B. The confinement force 

It is not easy to measure the free energy of confining 
the polymer in the nano-slit. Instead, the force that the 
confined polymer exerts on the plates is studied in the 
simulation. The results are compared with equation [3] 
from the theory. The simulation results for the force are 
plotted versus the theoretical relation in Fig. [51 A linear 
behavior is observed, which shows agreement between 


theory and simulation. 

The exponents 7 and w in the relation / oc N~^ D‘^ are 
obtained from simulation. According to Fig. [3 deviation 
in 7 from the theoretical value (= 1 ) is less than % 1 . 

In Fig. [5Ka), the value of u) is found from the logio- 
logio plot of force versus the nano-slit width. Deviation 
of the exponent from the theoretical value (-2.7) is less 
than % 8 . The simulation results are summarized in / oc 
]\[0-99jj-2.9^ Again, to check our definition for the nano- 
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log [(D+5)/b] 

FIG. 8. (Color online) logio-logio plot of / versus the nano¬ 
slit width, (a) and (b) show two different definitions for the 
nano-slit width (see Fig. n. As can be seen, omitting 5 results 
in an excellent improvement in the scaling exponent. 

slit width, the exponent uj is found from the logio-logio 
plot of force versus D + 6 (Fig. HKb)). Deviation of the 
exponent becomes more than %30. 


C. The effective spring constant 

Fluctuations in the radius of gyration (Ai?^) can be 
used to find an effective spring constant for the polymer, 
kef f ~ ■ This effective spring constant is a measure 

of the dynamics of the confined polymer. For example, 
the relaxation time of the confined polymer is found by 
dividing the friction coefficient from the solvent to this 
effective spring constant (T^ . 



FIG. 9. (Color online) Simulation results for the effective 
spring constant versus the theoretical relation N~^'^. O, 
□ , X, #, -I-, V, A, <1 and [> show the simulation results for 
D/b=l-9, in ascending order. 


The simulation results for the effective spring constant 
versus the theoretical relation (equation S]) are depicted 
in Fig. ini The linear curve confirms the excellent agree¬ 
ment between theory and simulation. This agreement is 
observed even when the nano-slit width is very small. 

The exponents f] and ( in the relation keff oc N^D'^ 
are found from the simulation results. Fig. [TOT a) shows 
the logiQ-logio plot of keff versus the polymer length. 
The average slope of the plots is -1.3. This exponent 
has the largest deviation, %13, among the exponents in¬ 
vestigated in this manuscript. The exponent is not a 
monotonic function of the nano-slit width and it cannot 
be deduced that an improved agreement between theory 
and simulation would be observed with increasing the 
nano-slit width. 

Dependence of the effective spring constant on the slit 
width is investigated in Fig. ITOT b). The average of the 
exponent is 0.48. Although deviation of the average from 
the theory is around %4, the exponents themselves are 
not very close to it. Overall, a good agreement between 
theory and simulation is observed and the simulation re¬ 
sults are summarized in the relation ke // (xiV-i-3D0-48. 

To check the effect of correlations in the calculation of 
(Ai?^|), a different averaging scheme was also examined. 
Data points with a time interval of IOOtq were collected 
from the whole data. Starting from different initial times, 
100 different series of data points can be defined. The 
values of the spring constant obtained from each of these 
series are very close to the ones shown on Figs. fTUf a) and 
IMb). 

Further investigation is needed to explain the larger de¬ 
viations of the dynamic exponents from theory, relative 
to the static ones. One can attribute these deviations to 
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log (D/b) 


FIG. 10. (Color online) (a) logio-logio plot of k^ff versus the 
polymer length, N, for different nanoslit width. O, □, x, •, 
+, V, A, <1 and l> show the simulation results for D/b=l-9, 
in ascending order. The corresponding rj values are -1.6, -1.6, 
-1.7, -1.4, -1.1, -1.5, -1.3, -1.3, -1.6. The average r; is -1.3. 
(b) logio-logio plot of kgff versus the slit width, for different 
polymer lengths. The average of ^ is 0.48. 


V. SUMMARY AND DISCUSSION 


In this manuscript, we checked the accuracy of the de 
Gennes’ blob model for a flexible polymer inside a nano¬ 
slit, using MD simulation. Theoretically, the blob model 
is accurate for wide slits and long polymers. However, 
we examined the model in the limit of very narrow slits. 
Nano-slits with widths equal to 1-9 times a monomer 
size were investigated. Simulation results both for the 
static properties and the dynamic spring constant of the 
polymer were in excellent agreement with the theory. 

Our results showed that the blob model can describe 
the effective spring constant of the confined polymer even 
for very narrow slits. This is in contrast to a recent 
study on single flexible polymers inside nano-channels. 
Indeed, the blob model is only applicable for the spring 
constant of the polymer in nano-channels with at least 
10 monomers width. Here, it was explained that this dis¬ 
crepancy is a result of the different origins of the polymer 
dynamics in the two cases. The internal dynamics of a 
blob is determining in the overall dynamics of the chain, 
for a polymer confined in a nano-channel. But, for a 
polymer confined in a nano-slit, the dynamics is domi¬ 
nated by the two-dimensional arrangement of the blobs 
between the plates. 

In our simulations, polymers with 100-300 monomers 
were used. It was observed that the effect of the finite 
size of the chains is not considerable. Our results on 
the accuracy of the blob model even for narrow slits and 
short chains can be very useful in the studies of polymers 
in confinement. Polymer confinement occurs in nano¬ 
technology devices, such as nano-pore sequencing 
DNA barcoding Q and polymer separation devices QT 

The main result of this manuscript was that the blob 
model is sometimes accurate beyond its theoretical limits, 
in practical situations. Considering that the blob model 
is used in many problems of polymer physics 2^ , this re¬ 
sult is applicable in different circumstances. A polymer 
adsorbed to a surface [1^ or compressed on a surface 


by an external field [2^ is confined to a region near the 
surface. For these problems with 2D confinement of the 
polymer, the above result can be used directly. How¬ 
ever, for a polymer under external tension more caution 
is needed, since the blobs are arranged in ID and the dy¬ 
namics may be more sensitive to the number of monomers 
inside one blob. 


the finite length of the chains or insufficient simulation 
times for equilibration of the dynamic quantities. How¬ 
ever, larger fluctuations may be an intrinsic property of 
the dynamic quantities and persist even in simulations 
with longer chains and longer run times. 
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